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ABSTRACT 


The  rotordynamic  response  of  an  imbalanced  rotor  accelerating  through  its  first 
lateral  bending  critical  speed  was  investigated  both  analytically  and  experimentally.  A  two 
degree-of-freedom  lumped  mass,  damping  and  stiffness  model  was  developed  to  simulate 
the  response  of  a  simply  supported,  single  disk  rotor  during  both  acceleration  and 
deceleration.  The  equations  of  motion  were  then  solved  numerically.  The  computer  model 
was  used  to  determine  the  effect  of  acceleration  rate,  asymmetric  stiffness  and  damping, 
and  acceleration  scheduling  on  the  maximum  amplitude  of  the  response.  Experimental  data 
for  a  simply  supported,  single  disk  rotor  accelerating  at  different  rates  were  compared  with 
the  computer  model.  Increased  acceleration  rates  and  damping  reduce  the  magnitude  of  the 
response.  Asymmetric  stiffness  and  acceleration  scheduling  can  also  be  used 
advantageously  to  reduce  the  maximum  amplitude  of  the  response. 
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I.  INTRODUCTION 


Rotating  machinery  is  used  widely  throughout  the  world  in  numerous  applications. 
These  applications  include  gas  and  steam  turbines  for  propulsion  and  electrical  power 
generation,  power  transmission  shafting,  compressors,  pumps,  gyroscopes,  high  speed 
flywheels,  and  several  other  devices  used  in  all  disciplines  and  areas  of  engineering. 

Many  of  these  rotating  machines  are  normally  operated  above  their  first  or  second 
critical  speeds  and  therefore  must  pass  through  one  or  more  critical  speeds  when  being 
started  and  again  when  coming  to  rest.  Therefore,  the  problem  of  determining  the 
maximum  amplitudes  and  stresses  in  the  system  due  to  the  forced  vibrations  during 
acceleration  through  critical  speed  is  of  vital  interest.  Of  even  greater  interest  is  the 
determination  of  what  methods  could  be  used  to  reduce  the  maximum  amplitudes  of  the 
lateral  bending  response  as  these  rotors  pass  through  their  critical  speeds.  This  is 
especially  important  in  applications  such  as  modem  turbomachinery  where  the  clearances 
between  the  rotor  and  the  housing  are  becoming  increasingly  tight.  The  current  trend 
toward  higher  spin  rates  to  raise  the  power  density  and  tighter  clearances  to  increase  the 
operating  efficiency  aggravates  the  problem  of  lateral  vibration  further.  As  the  demand  for 
high  technology  rotor  systems  increases,  so  does  the  need  for  the  ability  to  accurately 
predict  and  analyze  the  dynamic  response.  The  proper  prediction  and  understanding  of  the 
machine  dynamics,  especially  during  the  design  phase,  will  allow  for  the  use  of  lighter 
materials,  higher  obtainable  speeds,  and  more  power  dense  machines  while  decreasing 
maintenance  or  even  catastrophic  failure. 

Rotating  machinery  can  experience  several  modes  of  vibration.  The  lateral 
vibrations  are  of  primary  concern  here.  Due  primarily  to  residual  imbalance  in  the  rotating 
members,  lateral  vibrations  are  the  most  common  type  encountered.  Lateral  vibrations 
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bring  the  possibility  of  rubbing  or  equipment  damage  due  to  the  close  tolerances  of 
turbomachinery.  In  addition,  the  alternating  stresses  generated  by  these  vibrations  can 
damage  bearings,  and  non-synchronous  whirling  can  lead  to  fatigue  of  the  rotor  shafting. 

The  phenomena  associated  with  rotating  shafts  which  need  to  be  studied  and  fully 
understood  include;  the  response  at  or  transitioning  through  critical  speeds,  shaft  whirl,  oil 
whip,  chatter,  bearing  effects,  and  mounting  effects. 

The  objectives  of  this  research  were  to  investigate  how  the  different  physical 
parameters  such  as  rate  of  acceleration,  stiffness  and  damping  asymmetry,  and  the  effect  of 
acceleration  scheduling  effect  the  maximum  amplitude  of  the  response  as  a  rotor  is 
accelerated  through  its  first  lateral  bending  critical  speed.  In  addition,  to  further  develop 
the  rotordynamic  laboratory  at  the  Naval  Postgraduate  School  to  facilitate  the  study  of  an 
accelerating  rotor,  including  the  creation  of  a  mounting  table  for  isolation  of  the  rotor  and  to 
permit  mounting  of  the  rotor  kit  in  any  one  of  three  orientations  with  respect  to  gravity. 

In  Chapter  II  a  brief  background  of  the  pertinent  literature  on  non-constant 
frequency  excitation  to  support  the  research  is  provided. 

In  Chapter  III  a  two  degree-of-freedom  analytical  model  is  developed  in  which  the 
ability  to  vary  the  rate  of  acceleration,  amount  of  imbalance,  stiffness  and  damping  terms, 
and  initial  and  final  rpms  is  provided.  The  physical  characteristics  of  the  model  are 
discussed  and  the  equations  of  motion  are  provided  along  with  the  method  of  solution.  The 
results  of  the  model  are  compared  to  both  the  model  of  Lewis  [Ref.  1]  which  is  for  a  single 
degree-of-freedom  and  to  the  steady  state  Vejvoda  [Ref.  2]  model  given  an  extremely  slow 
rate  of  acceleration.  Parameters  are  varied  to  examine  how  each  effects  the  amplitude  of  the 
response  in  the  model  and  an  interpretation  of  the  results  is  provided. 

An  explanation  of  all  upgrades,  changes,  and  additions  to  the  experimental  facility 
is  provided  in  Chapter  IV  along  with  experimental  results  which  are  compared  to  the 
numerical  model. 
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Chapter  V  contains  a  discussion  of  the  results  from  both  the  numerical  models  and 


the  experiments. 

Finally,  conclusions  and  recommendations  are  provided  in  Chapter  VI. 
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II.  BACKGROUND 


A.  HISTORICAL  PERSPECTIVE 

The  first  documented  work  on  the  dynamics  of  rotating  shafts  was  presented  by 
Rankine  [Ref.  3]  in  1869.  His  analysis  dealt  with  the  steady  state  response  of  rotors. 
Neglecting  Coriolis  acceleration,  Rankine  came  to  the  mistaken  conclusion  that  supercritical 
operation,  that  is,  operating  at  an  angular  velocity  greater  than  one  or  more  of  the  systems 
natural  frequencies,  was  statically  unstable  and  therefore  not  possible  for  rotating 
machinery. 

In  1919,  Jeffcott  [Ref.  4]  published  what  has  become  one  of  the  definitive 
references  on  the  dynamics  of  rotating  shafts.  Analyzing  the  effects  of  imbalance  on  the 
steady  state  response  of  a  rotating  shaft,  Jeffcott  corrected  the  mistake  made  by  Rankine, 
showing  that  supercritical  operation  was  not  only  possible  but  could  be  beneficial. 
Obtaining  the  supercritical  condition,  however,  required  transitioning  through  the  first 
critical  speed  where  the  frequency  of  the  forcing  function  was  equal  to  the  first  lateral 
bending  natural  frequency  of  the  system  and  where  large  amplitude  vibrations  occurred. 
The  Jeffcott  rotor  was  developed,  which  is  comprised  of  a  long  shaft  with  a  concentrated 
mass  disk  and  provides  one  of  the  simplest  and  most  accurate  models  with  which  to  study  a 
simply  supported  rotating  shaft. 

The  first  analysis  which  dealt  with  a  system  accelerating  through  its  critical  speed 
was  published  in  1932  by  Lewis  [Ref.  1].  Lewis  found  the  exact  analytical  solution  of  a 
linearly  damped,  single  degree-of-freedom  system  accelerated  uniformly  from  rest  through 
its  critical  speed  in  terms  of  Fresnel’s  integrals.  The  figure  he  generated  for  the  case  of 
zero  damping  is  still  frequently  reproduced.  Lewis  verified  analytically  that  there  is  a 
decrease  in  the  amplitude  of  the  deflection  and  an  apparent  shift  in  the  critical  speed  of  a 
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system  which  accompanies  an  increase  in  the  rate  of  acceleration.  This  is  a  shift  to  the 
right,  or  a  higher  apparent  critical  speed,  for  acceleration  and  a  shift  to  the  left  for 
deceleration.  Lewis  also  attempted  to  provide  empirical  formulas  for  the  amount  of 
reduction  in  amplitude  and  the  amount  of  shift  given  large  values  of  acceleration.  Although 
his  model  was  not  one  of  a  rotating  reference  frame,  qualitatively  his  findings  are  still 
approximately  valid  for  a  rotating  shaft  experiencing  acceleration  through  its  critical  speed. 
The  Lewis  results  are  used  widely  as  rotors  are  normally  accelerated  through  their  critical 
speeds  as  quickly  as  possible. 

Study  in  the  area  of  acceleration  through  critical  speeds  continued.  In  1939,  Baker 
[Ref.  5]  published  his  results  on  the  Mathematical-Machine  determination  of  the  vibrations 
of  an  accelerated,  unbalanced  rotor. 

In  1947,  Meuser  and  Weibel  [Ref.  6]  presented  a  generalized  solution  for  the 
response  of  a  single  degree-of-freedom  system  having  linear  plus  cubic  elasticity,  and 
linear  or  viscous  damping  accelerating  through  its  critical  speed. 

In  1949,  McCann  and  Bennett  [Ref.  7]  extended  the  Lewis  model  to  a  two  degree- 
of-freedom  model  and  obtained  solutions  on  an  "electric-analog  computer." 

In  1966,  by  introducing  asymptotic  expressions  for  the  Fresnel  integrals,  Feam  and 
Millsaps  [Ref.  8]  derived  simple  algebraic  expressions  which  give  values  for  the  maximum 
amplitude  and  the  frequency  of  the  driving  force  where  the  maximum  displacement  occurs 
in  the  Lewis  model. 

Gluse  [Ref.  9]  in  1967,  investigated  the  case  of  a  rotor  being  accelerated  through  its 
critical  speed  by  a  constant  applied  torque  rather  than  a  constant  angular  acceleration.  His 
analysis  centered  on  the  case  of  an  unsuccessful  acceleration,  that  is,  an  acceleration  to  the 
vicinity  of  critical  speed  but  not  through  to  the  supercritical  area.  This  is  caused  by  the 
inertia  forces  which  produce  large  lateral  deflections  which  develop  counter  torques 
opposing  the  driving  torque  in  the  vicinity  of  the  system  natural  frequency.  Gluse 
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concluded  that  failure  to  accelerate  through  critical  speed  resulted  from  either  insufficient 
driving  torque  or  excessive  imbalance  of  the  rotor. 

In  1992,  Ishida,  Yamamoto,  and  Murakami  [Ref.  10]  investigated  analytically  the 
non-stationary  vibration  characteristics  during  acceleration  through  a  critical  speed  of  a  1/3- 
order  sub-harmonic  oscillation  of  forward  precession  of  a  rotating  shaft  system.  They 
concluded  that  the  maximum  amplitude  during  acceleration  depends  upon  the  angular 
acceleration  rate,  the  initial  disturbance  (if  one  exists),  the  initial  angular  velocity,  and  the 
initial  angular  position  of  the  imbalance. 

In  1988,  Vance  [Ref.  11]  analyzed  the  phenomena  caused  by  bearing  stiffness 
asymmetry  in  an  imbalanced  rotor.  His  observations  included  the  phenomena  of  a  split 
resonance  peak  in  the  response  caused  by  this  asymmetry  and  the  resultant  backward  whirl 
when  the  rotational  speed  is  between  these  two  frequencies. 

In  1994,  Vejvoda  [Ref.  2]  developed  a  two  degree-of-freedom,  steady  state 
analytical  model  and  utilized  an  experimental  rotor  facility  to  investigate  the  cause  of  split- 
resonance  and  backward  whirl  phenomena.  Isolating  the  effects  of  geometric  imperfections 
in  the  bearing  sleeves,  gravitational  forces,  and  bearing  support  stiffnesses,  Vejvoda 
concluded  that  the  stiffness  asymmetry  is  responsible  for  the  split  resonance  and  backward 
whirl  observed  in  his  experiments.  He  also  noted  that  the  rotor  system  orientation  with 
respect  to  the  gravitational  force  and  the  mount  isolation  caused  the  rotor  response  to 
change. 
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III.  ANALYTICAL  MODEL 


A.  MODEL  DEVELOPMENT 

In  order  to  predict  the  rotordynamic  response  of  a  rotor  accelerating  through  it's 
critical  speed,  an  analytical  model  was  developed.  The  rotor  was  modeled  as  a  two  degree- 
of-freedom  (x  and  y),  lumped  mass  and  lumped  stiffness  shaft  and  disk  assembly  as 
shown  in  Figure  1.  The  stiffnesses  and  cross-stiffnesses  of  the  shaft,  bearings,  and 
support  structures  have  been  combined  into  single  Kxx,  Kyy,  Kxy  and  Kyx  terms. 
Similarly,  the  proportional  damping  of  the  entire  system  has  been  combined  into  single 
Cxxj  Cyy,  Cxy,  and  Cyx  terms.  By  concentrating  all  of  the  systems  properties  into  a 
single  axial  plane,  this  model  provides  a  simple  way  to  parametrically  study  the  effects  of 
direct  and  cross-coupled  linear  stiffness  and  damping  on  the  maximum  response  amplitudes 
of  an  accelerating  rotor.  The  model  is  limited  in  that  no  variations  in  the  axial  direction  are 
taken  into  account,  only  forces  which  are  linearly  related  to  the  displacements  and  velocities 
of  the  rotor  are  included,  and  in  that  the  two  bearings  and  their  supports  of  the  simply 
supported  shaft  are  treated  as  identical,  but  not  necessarily  symmetrical  in  the  x  and  y 
directions.  Unlike  the  model  in  [Ref.  2]  in  which  only  the  steady  state  response  was  of 
interest,  an  acceleration  term  has  been  added  to  permit  the  investigation  of  rotor  response 
during  acceleration  and  deceleration. 

The  equations  of  motion  for  this  model  system  are  as  follows: 

M'x  +  C^x  +  C^y  +  K^x  +  K^y  =  mc(0cos 0  -  0^  sin 0)  (1) 

My  +  C^^y  +  Cy^x  +  K^^y  +  K^^x  =  me(0  sin  0+0^  cos  0)  (2) 
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Where  M  is  the  total  modal  mass  of  the  system  modeled  as  a  pinned-pinned  beam  as 
in  [Ref.  2],  m  is  the  mass  of  the  imbalance,  e  is  the  distance  from  the  center  of  the  shaft  to 
the  center  of  the  mass  imbalance  and  0is  constant.  In  the  steady  state  case,  the  angular 
acceleration  term  0  goes  to  zero  and  magnitude  of  the  forcing  function  is  dependent  only  on 
the  amount  of  imbalance  and  the  rotor  speed. 
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Figure  1.  Two  Degree-of-Freedom  Lumped  Mass  and  Stiffness  Model 
of  Rotating  Shaft  and  Disk  Assembly  with  Direct  and  Cross-Coupled  Linear 

Stiffness  and  Damping. 
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Equations  (1)  and  (2)  can  also  be  written  in  matrix  form  as: 


M  0 
0  M 


Cxx 

Cyx 


Kxy 

^yy 


•  2 

0  cos  0-0 
.  2 

0  sin  0  +  0 


sin  0 

COS0 


(3) 


B.  SOLUTION  TECHNIQUE 

These  equations  can  be  numerically  solved  using  the  Runge-Kutta  method.  For  this 
method  to  be  used,  the  two  second-order  differential  equations  (1)  and  (2)  can  be  written  as 
the  following  four  first-order  equations: 


>'1  =  ^2  (7) 


Several  computer  programs  were  written  in  MATLAB  to  solve  these  equations  with 
the  ability  to  enter  the  mass  imbalance,  Kxxj  ^yy»  ^xy>  Kyx>  Cxx>  ^yy*  ^xyj  Cyx,  initial 
and  final  shaft  RPM,  and  the  rate  of  acceleration  as  variables.  The  programs  utilized  the 
MATLAB  function  ODE45  which  uses  the  fourth  and  fifth  order  Runge-Kutta  methods. 
The  initial  conditions  for  this  numerical  method  were  determined  by  the  analytic  solution  of 
the  steady  state  (constant  speed)  response  for  the  initial  speed  selected.  9(t=0)  equal  to 
zero  was  used  but  could  be  changed  in  the  code. 

Since  the  ODE45  program  selects  it’s  own  time  divisions  for  the  numerical 
integration,  the  program  did  not  allow  for  the  interactive  entry  of  acceleration  as  a  function 
of  time  but  only  as  a  constant.  However,  acceleration  as  a  function  of  time  can  be  used  if 
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entered  directly  in  the  code  of  the  MATLAB  function  used  in  the  ODE45  program.  This 
function  contains  the  1st  order  differential  equations  of  motion  (equations  (4)  through  (7)) 
and  is  called  runga.m  for  the  acceleration  program.  In  addition,  a  program  was  written 
which  allowed  the  rate  of  acceleration  to  be  instantaneously  changed  at  a  given  RPM.  This 
allowed  for  the  limited  study  of  acceleration  scheduling  which  will  be  discussed  further  at 
the  end  of  this  chapter. 

C.  MODEL  VALIDATION 

In  order  to  validate  the  output  of  the  computer  models,  two  comparisons  were  made 
with  the  model  results: 

Lewis  [Ref.  1]  determined  the  exact  solution  to  the  problem  of  transitioning  a 
system  having  a  single  degree  of  freedom  and  linear  damping  through  its  critical  speed 
from  rest  at  a  uniform  acceleration.  The  plot  of  his  results  using  zero  damping  are  often 
reproduced.  The  Lewis  model  used  a  forcing  function  of  increasing  frequency  but  with  a 
constant  amplitude: 

f(t)  =  P  cos  (71  h  t2  +  c)  (8) 

where  P  is  the  magnitude  of  the  forcing  function,  h  is  the  acceleration  of  the  system  in 
cycles  per  second  per  second  and  a  is  the  initial  phase  angle  of  the  exciting  force. 

The  forcing  function  of  an  accelerating  imbalanced  rotor  is  not  only  one  of 
increasing  frequency  but  of  increasing  amplitude  as  well.  Therefore,  to  compare  the  model 
with  Lewis'  results,  it  was  changed  to  a  two  degree  of  freedom  system  with  forcing 
functions  in  the  x  and  y  directions  which  were  90  degrees  out  of  phase  and  were  of  a 
constant  amplitude. 

A  plot  of  the  results  showing  the  maximum  response  curves  for  various  values  of  q 
is  shown  in  Figure  2  in  the  non-dimensionalized  terms  used  by  Lewis.  The  term  r  is  the 
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number  of  undamped  free  vibrations  the  system  would  make  from  time  zero  to  time  t.  The 
non-dimensionalized  acceleration  rate  q,  is  given  by: 


?  = 


(9) 


where  N  is  the  natural  frequency  of  the  system  in  cycles  per  second  and  h  is  the 
acceleration  in  cycles  per  second  per  second.  The  term  q  can  also  be  stated  as  the  value  of  r 
at  which  the  instantaneous  frequency  of  the  exciting  force  is  equal  to  the  natural  frequency 
of  the  system.  The  amplitude  term  R,  is  given  by: 


R  = 


(10) 


where  x  is  the  displacement,  k  is  the  system  stiffness  and  P  is  the  magnitude  of  the  forcing 
function.  The  results  obtained  by  the  model  match  very  closely  with  the  analytical  solution 
presented  by  Lewis. 

The  next  comparison  made  was  with  the  steady  state  or  constant  excitation 
frequency  analytical  model  of  Vejvoda  [Ref.  2].  Since  the  ODE45  program  integrates  from 
an  initial  time  to  a  final  time,  zero  acceleration  could  not  be  entered  into  the  numerical 
program.  However,  an  extremely  slow  acceleration  should  fall  very  closely  to  the  results 
of  the  Vejvoda  model.  Figure  3  shows  the  amplitude  -vs-  rpm  plots  of  both  the  analytical 
steady  state  response  using  the  Vejvoda  model  and  the  numerical  model  response  with  an 
acceleration  of  1  rpm/sec.  This  plot  shows  that  the  numerical  model  with  a  very  slow  rate 
of  acceleration  does,  in  fact,  very  nearly  collapse  to  the  results  of  the  steady  state  analytical 
model. 
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XY  DISP  (mils) 


r/q  (non-dimensional  frequency) 


Figure  2.  Model  Results  Using  the  Lewis  Non-Dimensionalized  Terms 

with  Zero  Damping. 


SPEED  (rpm) 

Figure  3.  Comparison  of  the  Analytical  Model  with  an  Acceleration  Rate  of 
1  RPM/Sec  with  the  Vejvoda  Steady  State  Model. 


D.  MODEL  RESULTS 


1 .  Introduction 

Several  configurations  were  run  in  the  numerical  models  to  determine  the  effect  of 
each  of  the  various  parameters  on  the  maximum  amplitude  of  the  response.  In  all  cases 
presented,  an  imbalance  of  1  gram  at  a  distance  of  30.5  mm  from  shaft  centerline  was  used. 
The  system  mass  used  was  0.0061  Ibf-sec^/in,  which  is  the  modal  mass  determined  for  the 
experimental  facility.  When  stiffness  coefficients  were  kept  constant  they  were  in  the 
region  of  495  to  515  Ibf/in  which  is  the  approximate  modal  stiffness  values  determined  for 
the  experimental  facility.  Likewise,  when  damping  was  kept  constant  it  was  held  at  a  value 
around  0.04  Ibf-sec/in  which  gives  a  realistic  value  of  the  damping  ratio  ^  which  matches 
experimental  data.  Each  of  the  parameters  was  varied  individually  to  determine  its  effect  on 
the  response.  The  goal  was  to  determine  what  factors  reduced  the  maximum  response 
amplitude  and/or  the  duration  or  number  of  cycles  experiencing  high  amplitudes. 

The  displacement  curves  shown  in  the  plots  are  the  maximum  radial  displacement, 
XY  given  by  the  formula: 


xr 


(11) 


The  term  XY  is  used  in  lieu  of  the  more  common  terms  for  radial  displacement,  r  or 
R  to  avoid  confusion  with  the  Lewis  non-dimensional  terms  (see  Figure  4). 
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X  DISP  (mils) 


Figure  4.  Diagram  Showing  the  Relationship  between  the  Radial 
Displacement  XY  and  the  X  and  Y  Displacement. 

2.  Determination  of  number  of  cycles  experienced  between 
0.9(D/a)n  and  l.lto/ojn. 

One  of  the  goals  was  to  minimize  the  number  of  cycles  experienced  at  the  high 
vibration  levels.  To  determine  the  number  of  cycles  occurring  between  plus  or  minus  10% 
of  ®/Q)n  =  1,  the  following  formula  was  derived: 


Since, 


and  the  range  of  to  is  from  0.9  to  1.1  ton  , 


Then, 


2n 


(in  hertz) 


(12) 


(13) 
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and  the  time  required  to  traverse  this  range  is  given  by: 


h  2%h 


(14) 


Where  h  is  the  rate  of  acceleration  in  cycles  per  second  per  second. 

The  number  of  cycles  experienced  while  traversing  this  range  can  be  expressed  as: 


1  2 

#  Cycles  =  0)  t  +  —ht 


(15) 


Using  substitution,  the  resultant  equation  is: 


#  Cycles  = 


0.05  K 
n^Mh 


(16) 


To  put  this  in  terms  of  q,  use  the  relationships: 


(17) 


where  N  is  the  natural  frequency  in  cycles/sec,  and: 


This  gives  us  a  very  simple  form  of  the  equation: 


#  Cycles  =  0.2q  (19) 

Equations  (13)  and  (16)  give  us  the  number  of  cycles  that  occur  within  plus  or 
minus  10%  of  to/(On  =  1.  At  high  rates  of  acceleration,  the  shift  in  the  apparent  critical 
speed  can  place  the  region  of  highest  amplitudes  out  of  this  range.  The  number  of  cycles 
occurring  within  the  same  range  around  the  ®  of  maximum  response  will,  however,  be 
very  close  to  this  same  value. 

3.  Symmetric  stiffness  and  damping  at  different  rates  of 
acceleration. 

The  maximum  displacement-vs-rpm  curves  for  various  acceleration  rates  on  a 
symmetric  rotor  are  shown  in  Figure  5.  Kxx  and  Kyy  are  symmetric  at  495  Ibf/in  and  Cxx 
and  Cyy  are  equal  at  0.04  Ibf-sec/in  (^=.0115).  All  cross  terms  are  zero.  Acceleration 
rates  of  1,  250,  500,  1000,  and  2000  rpm/sec  were  used  (these  acceleration  rates  have 
approximate  non-dimensional  q  values  of  123200,  493,  246,  123,  and  62  respectively). 
Just  as  in  the  Lewis'  curves,  these  curves  show  the  maximum  amplitude  decreasing  with 
increasing  acceleration  rates.  The  apparent  shift  of  the  critical  speed  is  also  easily  seen  in 
this  figure.  As  the  acceleration  is  doubled  from  250  rpm/sec  to  500  rpm/sec  there  is  a 
12.9%  decrease  in  the  maximum  amplitude  and  a  40  rpm  increase  in  the  speed  of  maximum 
response.  As  acceleration  is  again  doubled  to  1000  rpm/sec,  there  is  an  accompanying 
15.5%  decrease  in  the  maximum  amplitude  and  a  increase  of  60  rpm.  When  the 
acceleration  rate  is  again  doubled  to  1000  rpm/sec,  the  maximum  amplitude  decreases  by 
19.5%  and  the  increase  in  the  speed  of  maximum  response  is  80  rpm. 
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4.  Symmetric  stiffness  and  damping  at  different  rates  of 
deceleration. 

The  same  values  for  all  parameters  used  in  the  accelerating  case  were  again  used  for 
decelerating  at  the  same  rates.  Figure  6  shows  the  corresponding  curves  for  the  same  rates 
of  deceleration.  This  shows  slightly  lower  maximums  than  for  the  accelerating  case.  The 
maximum  amplitude  is  0.03%  lower  for  the  1  rpm/sec  curve  and  increases  to  6.4%  lower 
for  the  1000  rpm/sec  curve.  The  trend  then  reverses  and  the  curve  for  2000  rpm/sec 
deceleration  is  only  2.9%  less  then  its  corresponding  acceleration  curve.  The  shift  in  the 
apparent  critical  speed  is  now  to  the  left  and  by  an  amount  similar  to  that  seen  in  the  case  of 
acceleration. 

5.  Effects  of  stiffness  asymmetry 

To  determine  if  stiffness  asymmetry  could  be  used  advantageously  to  reduce  the 
maximum  amplitude,  model  rans  with  varying  degrees  of  stiffness  asymmetry  were  made. 
The  effect  of  varying  the  stiffness  asymmetry  on  a  rotor  experiencing  acceleration  is  shown 
in  Figure  7.  The  acceleration  for  all  curves  was  held  constant  at  360  rpm/sec  (q=356),  Cxx 
and  Cyy  were  held  at  0.04  Ibf-sec/in  (^=.0113),  and  the  stiffness  asymmetry  was  varied 
from  a  delta  of  zero  to  70  Ibf/in  in  steps  of  10  Ibf/in  around  a  mean  of  5 15  Ibf/in.  All  cross 
terms  were  zero.  With  the  stiffness  symmetric  at  515  Ibf/in,  the  maximum  amplitude  is 
around  16  mils.  As  a  slight  asymmetry  of  10  Ibf/in  is  introduced,  the  maximum  jumps  to 
around  19  mils,  an  increase  of  almost  19%.  As  the  asymmetry  is  increased  to  20  Ibf/in,  the 
maximum  increases  to  almost  21  mils  or  by  31%.  As  the  asymmetry  is  further  increased, 
the  maximum  now  starts  to  decrease  due  to  separation  of  the  split  resonance  peaks.  This 
split  resonance  was  determined  by  Vejvoda  [Ref.  2]  to  be  one  cause  of  backwhirl.  At  an 
asymmetry  of  70  Ibf/in  the  maximum  amplitude  has  returned  very  nearly  to  the  value  at  zero 
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asymmetry,  however,  there  are  now  two  distinct  critical  speeds  and  a  wider  range  of  rpms 
in  which  the  amplitude  of  the  response  is  significant. 

To  determine  if  the  trend  is  the  same  in  the  steady  state  case  of  zero  acceleration,  the 
same  values  were  run  in  the  Vejvoda  model  and  the  results  are  shown  in  Figure  8.  The 
trend  appears  to  be  the  same  regardless  of  whether  the  rotor  is  accelerating  or  not. 

6.  Effects  of  a  cross  stiffness  term. 

In  order  to  determine  the  effects  of  a  cross  stiffness  term,  acceleration  was  held 
constant  at  360  rpm/sec  (q=356),  Cxx  and  Cyy  at  0.04  Ibf-sec/in  (^=.01 13),  Kxx  and  Kyy 
were  held  constant  at  515  Ibf/in,  and  Kyx  was  varied  from  zero  to  30  Ibf/in.  Kxy  and  the 
damping  cross  terms  were  held  at  zero.  The  results  are  shown  in  Figure  9.  The  cross 
stiffness  term  appears  to  have  no  significant  effect  on  the  maximum  amplitude  of  the 
accelerating  rotor  until  it  reaches  a  value  of  30  Ibf/in.  It  then  has  a  very  significant  effect, 
increasing  the  maximum  by  nearly  40%  as  well  as  shifting  the  apparent  critical  speed  to  the 
right.  In  addition,  starting  at  a  value  of  Kyx  =  10  Ibf/in,  this  cross  term  induced  backwhirl 
in  the  accelerating  rotor  model.  The  rotor  orbit  from  2860  to  2875  rpm  showing  the  onset 
of  backwhirl  is  shown  in  Figure  10. 

Again,  to  determine  if  the  trend  is  the  same  in  the  psuedo-static  case  of  zero 
acceleration,  the  same  values  were  ran  in  the  Vejvoda  model  and  the  results  are  shown  in 
Figure  11.  In  the  case  of  zero  acceleration,  the  results  were  very  different.  Unlike  in  the 
accelerating  case,  the  cross  term  introduced  in  the  pseudo-static  case  had  a  significant  effect 
on  the  maximum  amplitude  for  each  value  of  Kyx-  As  the  value  of  the  cross  term 
increased,  so  did  the  amplitude  of  the  response.  However,  this  did  not  induce  backwhirl  in 
the  static  model  and  produced  no  significant  shift  in  the  apparent  critical  speed. 
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An  additional  run  was  made  was  made  with  an  acceleration  rate  of  16.67  rpm/sec 
(q=7689)  and  as  can  be  seen  in  Figure  12,  these  curves  begin  to  show  the  same  trend  as  in 
the  case  of  zero  acceleration. 

7.  Effects  of  symmetric  damping. 

Although  its  effects  are  well  documented,  symmetric  damping  was  evaluated  for  an 
accelerating  rotor  in  the  numeric  model  and  the  results  are  shown  in  Figure  13.  The 
acceleration  rate  was  held  constant  at  360  rpm/sec  (q=342),  Kxx  and  Kyy  were  held  at  495 
Ibf/in  and  symmetric  damping  was  evaluated  at  zero,  0.02,  0.04,  and  0.08  (^=0,  0.0058, 
0.01 15  and  0.0230).  All  cross  terms  were  zero.  As  expected  from  the  Lewis  analysis,  the 
amplitude  of  the  maximum  response  decreases  significantly  with  an  increase  in  damping. 

8.  Effects  of  damping  asymmetry. 

The  effect  of  varying  the  degree  of  damping  asymmetry  on  a  rotor  experiencing 
acceleration  is  shown  if  Figure  14.  The  acceleration  for  all  curves  was  held  constant  at  360 
rpm/sec  (q=356),  Kxx  and  Kyy  were  held  at  515  Ibf/in,  and  the  damping  asymmetry  was 
varied  from  a  delta  of  zero  to  0.05  Ibf-sec/in  in  steps  of  0.01  Ibf-sec/in  around  a  mean  of 
0.04  Ibf-sec/in  (^=0.01 13).  All  cross  terms  were  zero.  With  each  increase  in  the  amount 
of  damping  asymmetry,  there  was  an  accompanying  increase  in  the  amplitude  of  the 
maximum  response  of  6  to  10%.  This  maximum  value  was  the  same  as  what  would  be 
obtained  by  having  symmetric  damping  at  the  lower  value.  That  is,  if  Cxx  =  0.065  and 
Cyy  =  0.015,  the  resulting  maximum  amplitude  is  the  same  as  would  be  observed  with 
Cxx  =  Cyy  =  0.015.  This  trend  did  not  reverse  as  it  did  in  the  asymmetric  stiffness  case. 
This  asymmetric  damping  also  induced  backwhirl.  The  rotor  orbit  from  2885  to  2900  rpm 
displaying  the  onset  of  the  backwhirl  is  shown  in  Figure  15. 

When  compared  with  the  zero  acceleration  runs  made  in  the  Vejvoda  model  (see 
Figure  16),  the  trend  appears  to  be  the  same.  However,  there  is  no  induced  backwhirl  in 
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the  zero  acceleration  model.  It  therefore  appears  that  damping  asymmetry  can  induce 
backwhirl  in  an  accelerating  rotor  but  not  in  a  steady  state  rotor. 

9.  Effects  of  acceleration  changes. 

While  studying  the  responses  from  the  numerical  model  at  different  rates  of 
acceleration,  the  question  arose  as  to  whether  acceleration  scheduling  could  be  used  to 
minimize  the  response  of  the  rotor.  The  possibility  of  taking  advantage  of  the  apparent 
shift  in  the  critical  speed  at  high  rates  of  acceleration  was  explored.  The  idea  was  to  start  at 
a  high  rate  of  acceleration  which  would  cause  a  shift  in  the  critical  speed  to  the  right.  After 
the  rotor  is  past  what  would  be  the  critical  speed  at  a  lower  rate  of  acceleration,  but  before 
the  rotor  reaches  the  critical  speed  at  its  present  higher  acceleration,  the  acceleration  would 
be  immediately  slowed  to  the  lower  rate.  The  question  was,  would  this  "fool"  the  rotor 
into  thinking  that  it  was  already  past  its  critical  speed  and  therefore  would  the  amplitude 
start  decreasing  to  give  a  maximum  less  than  what  it  would  have  been  if  the  acceleration 
were  left  at  the  higher  rate.  Figure  17  shows  this  more  clearly. 

A  numerical  model  was  modified  to  allow  for  one  instantaneous  change  in 
acceleration  with  the  rpm  at  which  the  change  was  to  occur  entered  as  a  variable.  This 
model  was  limited  in  that  it  gave  an  unrealistic  instantaneous  change  in  acceleration,  that  is: 

e  =  9^(t<to)  (20) 

e  =  e^{t>tj  (21) 


and  that  it  did  not  allow  for  any  complex  scheduling  of  acceleration  in  that  0(t)  was  not 
specified. 

Figure  18  shows  the  results  of  one  of  these  change  acceleration  runs.  The  stiffness 
values  were  Kxx  and  Kyy  =  495  and  the  damping  values  were  Cxx  and  Cyy  =  0.04 
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(^=0.01 15).  The  initial  acceleration  was  2000  rpm/sec  (q=62)  which  was  switched  to  250 
rpm/sec  (q=493)  at  an  rpm  of  2870.  The  figure  shows  the  results  of  the  run  along  with 
plots  of  the  250  and  2000  rpm/sec  acceleration  response  curves.  As  can  be  seen,  the 
amplitude  of  the  response  is  only  slightly  less  than  that  of  the  continuous  2000  rpm/sec 
acceleration  response.  However,  the  response  does  "ring  down"  much  more  rapidly 
yielding  fewer  cycles  at  the  higher  vibration  levels. 
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Figure  5.  Symmetric  Stiffness  and  Damping  at  Various  Acceleration  Rates 
(Kxx=Kyy=495  Ibf/in,  Cxx=Cyy=0.04  Ibf-sec/in,  All  Cross  Terms  are 

Zero). 
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Figure  6.  Symmetric  Stiffness  and  Damping  at  Various  Deceleration  Rates 
(Kxx=Kyy=495  Ibf/in,  Cxx=Cyy=0.04  Ibf-sec/in,  All  Cross  Terms  are 

Zero). 
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Figure  7.  Effects  of  Stiffness  Asymmetry  around  a  Mean  of 
515  Ibf/in  (Cxx=Cyy=0.04  lbf>sec/in,  All  Cross  Terms  are  Zero, 
Acceleration  Rate  is  360  RPM/Sec). 
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Figure  8.  Effects  of  Stiffness  Asymmetry  around  a  Mean  of 
515  IbMn  in  the  Steady  State  Model  (Cxx=Cyy=0.04  Ibf-sec/in,  All  Cross 

Terms  are  Zero). 
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Figure  9*  Effects  of  a  Stiffness  Cross  Term  (Cxx=Cyy=0.04  Ibf-sec/in, 
Kxy  and  Damping  Cross  Terms  are  Zero,  Acceleration  Rate  is  360 

RPM/Sec). 


Figure  10.  Rotor  Orbit  Showing  the  Onset  of  Backwhirl  between  2860  and 
2875  RPM  due  to  a  Stiffness  Cross  Term  (Cxx=Cyy=0.04  Ibf-sec/in,  Kxy 
and  Damping  Cross  Terms  are  Zero,  Acceleration  Rate  is  360  RPM/Sec). 
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Figure  11.  Effects  of  a  Stiffness  Cross  Term  in  the  Steady  State  Model 
(Cxx=Cyy=0.04  Ibf-sec/in,  Kxy  and  Damping  Cross  Terms  are  Zero). 
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Figure  12.  Effects  of  a  Stiffness  Cross  Term  (Cxx=Cyy=0.04  Ibf-sec/in, 
Kxy  and  Damping  Cross  Terms  are  Zero,  Acceleration  Rate  is  16.67 

RPM/Sec). 
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Figure  13.  Effects  of  Symmetric  Damping  (Kxx=Kyy=495  Ibf/in,  all  Cross 
Terms  are  Zero,  Acceleration  Rate  is  360  RPM/Sec). 
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Figure  14.  Effects  of  Damping  Asymmetry  around  a  Mean  of  0.04  Ibf- 
sec/in  (Kxx=Kyy=515  Ibf/in,  All  Cross  Terms  are  Zero,  Acceleration  Rate 

is  360  RPM/Sec). 


X  DISP  (mils) 


Figure  15.  Rotor  Orbit  Showing  the  Onset  of  Backwhirl  between  2885  and 
2900  RPM  due  to  Damping  Asymmetry  (Kxx=Kyy=515  Ibf/in,  All  Cross 
Terms  are  Zero,  Acceleration  Rate  is  360  RPM/Sec). 
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Figure  16.  Effects  of  Damping  Asymmetry  around  a  Mean  of  0.04  Ibf- 
sec/in  in  the  Steady  State  Model  (Kxx=Kyy=515  Ibf/in,  All  Cross  Terms 
are  Zero,  Acceleration  Rate  is  360  RPM/Sec). 
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Figure  17.  Representation  of  how  the  Shift  in  the  Apparent  Critical  Speed 
may  he  used  to  reduce  Maximum  Rotor  Response. 


Figure  18.  Results  of  Instantaneously  Changing  the  Acceleration  Rate 
(Kxx=Kyy=495  Ihf/in,  Cxx=Cyy=0.04  Ihf-sec/in,  All  Cross  Terms  are 

Zero). 
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IV.  EXPERIMENTAL  INVESTIGATION 


A.  EXPERIMENTAL  FACILITY 

1.  Introduction 

The  experimental  facility  used  consisted  of  the  Bentley  Nevada  rotor  kit  and 
Labview  software  used  by  both  Simei  [Ref.  12]  and  Vejvoda  (Ref.  2].  However,  several 
updates  and  or  changes  were  made  to  the  facility.  Figures  19  through  21  show 
photographs  of  the  facility. 

2.  Upgrade  to  Pentium 

The  DataStor  IBM  personal  computer  used  for  data  acquisition  was  replaced  with  a 
Unitek  personal  computer  with  a  Pentium  chip  and  a  CPU  speed  of  66MHz.  The  National 
Instmments  Data  Acquisition  Interface  Board  (AT-MIO-16F-5)  is  the  same. 

3.  Table  Mounting 

In  both  Simei's  and  Vejvoda's  experiments,  the  rotor  kit  was  placed  on  a  foam  pad 
sitting  on  an  ordinary  table.  In  Vejvoda’s  experiments  with  the  kit  mounted  in  the  vertical 
and  wall  mounted  horizontal  planes,  the  rotor  kit  was  mounted  to  a  non-structural  wall  of 
the  laboratory.  The  vibrational  characteristics  of  these  mountings  and  therefore  the  effects 
they  had  on  the  rotor  response  were  unknown.  It  is  desirable  to  have  the  capability  to 
firmly  mount  the  rotor  kit  in  any  one  of  the  three  orientations  with  respect  to  gravity  on  a 
table  with  a  natural  frequency  which  is  much  lower  than  the  frequency  of  the  forcing 
function  which  might  be  experienced  by  the  imbalanced  rotor.  In  this  way,  the  mounting 
table  will  not  be  excited  by  the  imbalanced  rotor  and  therefore,  will  not  effect  the  response 
of  the  rotor. 
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A  heavy  steel  table  was  designed  for  the  rotor  lab  and  is  shown  in  Figures  19, 20, 
and  21  with  the  rotor  kit  placed  in  the  flat  horizontal,  wall  mounted  horizontal,  and  vertical 
mounting  positions  respectively.  The  weight  of  the  table  is  266  lbs  without  the  rotor  kit 
and  it  is  supported  by  four  Barry  6E100  machinery  mounts  which  have  a  very  low  stiffness 
(each  rated  at  5.4  to  8.1  Hz  for  loads  of  50  to  100  Ibf).  This  provides  for  a  natural 
frequency  of  less  than  8  Hz  in  both  the  axial  and  radial  directions. 

4.  Optical  Encoder  and  Frequency  Divider 

The  Labview  program  ROTOR. VI  used  by  Simei  and  Vejvoda  used  a  data 
acquisition  scan  rate  of  4096  scans/sec.  This  scan  rate  was  constant  regardless  of  rotor 
speed,  therefore  the  exact  position  of  the  shaft  when  each  acquisition  was  made  was  not 
controlled.  It  is  preferable  to  know  the  position  of  the  shaft  at  each  acquisition  for  the 
purpose  of  phase  locking  and  ensemble  averaging.  To  accomplish  this,  a  Hewlett  Packard 
HEDS5540  optical  encoder  was  attached  to  the  end  of  the  drive  motor  shaft.  This  encoder 
outputs  256  pulses  per  shaft  revolution.  This  output  is  sent  to  a  frequency  divider  which 
was  installed  on  the  National  Instrument  general-purpose  termination  breadboard  (SC- 
2070).  This  frequency  divider  takes  the  256  pulse/rev  signal  and  converts  it  to  either  a 
128,  64,  32,  or  16  pulse/rev  signal.  This  output  is  then  connected  to  the  breadboard's 
OUT2  terminal  to  be  used  as  an  external  clocking  signal.  The  ROTOR. VI  program  has 
been  modified  to  use  this  external  clocking  signal  for  data  acquisition  vice  the  internal  scan 
rate.  A  complete  wiring  diagram  of  the  experimental  facility  which  includes  the  frequency 
divider  connections  for  the  desired  scans  per  revolution  is  shown  in  Figure  22. 

5.  Additional  Labview  Programs 

Some  additions  and  modifications  were  made  to  the  Labview  data  acquisition 
programs  to  enhance  the  usability  of  the  lab. 
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ROTOR2.VI  is  a  program  identical  to  the  ROTOR. VI  described  in  [Ref.  12]  with 
the  exception  of  writing  the  data  to  a  file.  Placing  this  program  into  continuous  run  allows 
the  user  to  observe  real  time  data  on  the  running  rotor  kit  without  filling  up  a  data  file.  This 
can  be  used  much  like  an  oscilloscope  to  observe  the  current  rotor  orbit  and  y  displacement 
versus  time  plots. 

ACC  ROTOR. VI  is  a  Labview  program  which  allows  data  acquisition  on  an 
accelerating  rotor.  This  program  basically  places  the  ROTOR. VI  program  into  a  while  loop 
but  has  also  added  a  stop  button  and  a  filename  input  to  the  control  panel.  This  program 
will  continue  to  write  data  to  a  file  during  acceleration  until  the  stop  control  is  activated. 
This  is  the  program  that  was  used  to  take  the  experimental  data  discussed  later  in  this 
chapter. 

RPM.VI  is  a  Labview  program  which  utilizes  the  AM9513  clock  circuits  in  the  data 
acquisition  board  to  determine  the  frequency  of  the  external  clocking  signal.  Since  the 
optical  encoder  provides  a  set  number  of  clocking  signals  per  revolution,  the  program  can 
determine  the  shaft  rpm  and  thereby  eliminate  the  need  for  the  stroboscope.  This  is  very 
useful  in  instantaneously  knowing  the  approximate  shaft  speed  and  acceleration  but  it  was 
found  that  the  last  digit  of  the  readout  changes  at  such  a  rapid  rate  that  it  is  difficult  to  read 
the  exact  rpm. 

COUNTER  SUB  .VI  is  a  Labview  program  which  similarly  uses  the  AM9513  clock 
circuits  to  determine  the  time  interval  between  the  external  clocking  signals  and  write  the 
interval  to  a  file.  This  program  was  written  to  be  used  in  conjunction  with  the  ACC 
ROTOR. VI.  It  was  desired  to  be  able  to  obtain  x,  y,  and  time  interval  data  during  an 
accelerating  run.  In  this  way,  the  actual  acceleration  rate  could  be  determined  and  the  exact 
rpm  at  each  acquisition  could  be  determined.  Attempts  to  merge  these  programs  or  to  run 
them  simultaneously  were  unsuccessful  using  Labview  version  3.1,  probably  due  to  both 
programs  trying  to  access  the  same  registers  simultaneously. 
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Figure  19.  Photograph  of  Lab  Facility  with  the  Rotor  Mounted  in  the  Flat 

Horizontal  Position. 
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Figure  20.  Photograph  Shoiving  the  Rotor  Kit  in  the  Wall  Mounted 

Horizontal  Position. 


Figure  21.  Photograph  Showing  the  Rotor  Kit  Mounted  in  the  Vertical 

Position. 
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Figure  22.  Experimental  Facility  Wiring  Diagram  and  Frequency  Divider 

Connections. 


B.  EXPERIMENTAL  RESULTS 


1.  Analysis  of  the  table  mounting 

The  first  experiments  conducted  with  the  rotor  kit  mounted  on  the  new  table 
consisted  of  rerunning  some  of  those  made  in  [Ref.  2],  Vejvoda  had  determined 
approximate  stiffness  and  damping  coefficients  for  all  three  planes  of  reference 
experimentally.  To  determine  if  his  mounting  had  affected  those  values,  or  if  the  present 
table  mounting  changed  the  values,  the  zero  acceleration  runs  in  all  three  planes  were 
repeated.  The  experimental  configuration  was  the  same  as  shown  in  Figure  5  of  [Ref.  2]. 
Figure  23  shows  the  results  of  the  flat  horizontal  run  using  1  gram  imbalance  and  a  plot  of 
the  Vejvoda  analytical  model  results  for  the  K  and  C  values  given.  Figures  24  and  25 
show  the  same  for  the  wall  mounted  and  vertical  runs  respectively.  The  stiffness  and 
damping  coefficients  have  changed  considerably  with  the  rotor  kit  mounted  on  the  heavy, 
resiliently  mounted  table.  In  addition,  there  was  no  backwhirl  observed  in  the  flat 
horizontal  position  as  was  observed  in  Vejvoda's  experiments. 

2.  Constant  Acceleration  runs 

The  experimental  set-up  was  not  as  flexible  as  the  numerical  model  in  that  the 
stiffness  and  damping  values  could  not  be  readily  controlled.  However,  the  motor 
controller  for  the  rotor  kit  did  allow  for  control  of  constant  acceleration  rates  up  to  15000 
rpm/min  within  an  advertised  accuracy  of  10  percent.  Several  runs  were  made  in  the  flat 
horizontal  mounting  position  using  1  gram  imbalance  and  acceleration  rates  of  3000, 6000, 
9000,  12000,  and  15000  rpm/min.  Figures  26  through  30  show  experimental  results 
against  the  numerical  predictions  using  the  K  and  C  values  determined  in  the  previous 
paragraph.  Since  the  rpm  at  which  the  response  occurred  in  the  experimental  data  could  not 
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be  determined  using  our  present  data  acquisition  system,  the  maximum  amplitudes  are  what 
is  of  interest.  Some  results  are  tabulated  in  Table  1. 


Table  1.  Sample  of  Experimental  Results 


Acceleration 

Rate  (rpm/min) 

Max  Amplitude 

of  Model  (mils) 

Max  Experimental 

Amplitude  (mils) 

%  Error 

3000 

13.29 

13.41 

<  1% 

6000 

13.21 

13.3 

<  1% 

9000 

13.05 

13.09 

<  1% 

12000 

12.87 

13.09 

1.7% 

15000 

12.67 

12.77 

<  1% 

The  experimental  values  of  the  maximum  amplitude  obtained  were  slightly  higher 
than  in  the  model  for  all  rates  of  acceleration.  The  decreasing  trend  in  the  amplitude  of  the 
maximum  response  as  acceleration  is  increased  matches  very  closely  with  the  numerical 
model.  The  experimental  data  was  within  2%  of  the  analytical  in  all  cases  and  varied  less 
than  1%  for  most. 


40 


XY  DISP  (mils)  Y  DISP  (mils)  .  X  DISP  (mils) 


2200  2400  2600  2800  3000  3200  3400 


2200  2400  2600  2800  3000  3200  3400 

SPEED  (rpm) 

20 1  I  I  I - 1 - , - 1 


0^  - 1 - 1 - 1 _ I _ t  -1 

2200  2400  2600  2800  3000  3200  3400 


SPEED  (rpm) 


Figure  23.  Experimental  and  Theoretical  X  Displacement,  Y  Displacement, 
and  XY  Displacement  Steady  State  Data  (Kxx=493,  Kyy=510,  Cxx=0.07, 
Cyy=0.088,  All  Cross  Terms  are  Zero  (Flat  Horizontal  Mounting)). 
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X  Displacement,  Y  Displacement, 
(Kxx=493,  Kyy=497,  Cxx=0.08, 
(Wall  Mounted  Horizontal)). 
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Figure  25.  Experimental  and  Theoretical  X  Displacement,  Y  Displacement, 
and  XY  Displacement  Steady  State  Data  (Kxx=521,  Kyy=503,  Cxx=0.068, 
Cyy=0.084,  All  Cross  Terms  are  Zero  (Vertical  Mounting)). 
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Figure  26.  Experimental  -vs-  Theoretical  Maximum  XY  Displacement  Plot 
(Acceleration  rate  is  15000  RPM/Min). 


RPM 


Figure  27.  Experimental  -vs-  Theoretical  Maximum  XY  Displacement  Plot 
(Acceleration  rate  is  12000  RPM/Min). 
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Figure  28.  Experimental  -vs-  Theoretical  Maximum  XY  Displacement  Plot 
(Acceleration  rate  is  9000  RPM/Min). 
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Figure  29.  Experimental  -vs-  Theoretical  Maximum  XY  Displacement  Plot 
(Acceleration  rate  is  6000  RPM/Min). 
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Figure  30.  Experimental  -vs-  Theoretical  Maximum  XY  Displacement  Plot 
(Acceleration  rate  is  3000  RPM/Min). 
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V.  DISCUSSION  OF  RESULTS 


Both  the  computer  model  and  the  experimental  data  reveal  that  increases  in  the  rate 
of  acceleration  and  the  amount  of  damping  have  the  greatest  effect  on  reducing  the 
amplitude  of  the  response  of  the  rotor  passing  through  it  first  lateral  bending  critical  speed. 
The  rate  of  acceleration  must  be  relatively  large  to  have  a  significant  effect  on  reducing  the 
maximum  amplitude  of  the  response.  Lewis  [Ref.  1]  considered  a  q  value  of  500  to  still  be 
a  rather  slow  rate  of  acceleration.  This  analysis  showed  in  Figure  5  that  a  q  value  of  almost 
500  still  gave  a  noticeable  reduction  of  greater  than  12  percent  in  the  amplitude  of  the 
maximum  response.  Therefore  the  acceleration  should  give  a  q  value  of  less  than  500  to 
have  notable  effects.  If  these  large  accelerations  can  be  achieved,  they  are  possibly  the 
easiest  way  to  significantly  reduce  the  magnitude  of  the  response. 

Acceleration  introduces  other  phenomena  as  well.  In  addition  to  the  apparent  shift 
in  the  critical  speed  during  acceleration,  the  analysis  also  uncovered  that  acceleration  can 
generate  backwhirl  under  conditions  in  that  it  would  not  be  generated  with  little  or  no 
acceleration.  This  was  predicted  by  the  model  in  both  the  case  of  a  cross-stiffness  term  and 
in  the  case  of  direct  damping  asymmetry.  Another  effect  of  acceleration  is  that  the  whirling 
is  not  synchronous  and  therefore  the  phase  shift  of  the  mass  imbalance  transitioning  from 
the  outside  to  the  inside  of  the  orbit  as  the  rotor  passes  through  critical  is  no  longer  evident. 
This  non-synchronous  whirling  can  lead  to  fatigue  in  the  rotor  shaft.  The  benefit  of 
reducing  the  maximum  response  due  to  a  high  rate  of  acceleration,  however,  outweighs  the 
detriments. 

Increased  damping  is  another  highly  effective  means  to  keep  the  relative  maximum 
amplitudes  low,  though  it  is  much  more  difficult  to  predict  and  to  control.  The  assessment 
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of  asymmetric  damping  revealed  that  for  lateral  bending,  the  total  damping  is  only  as  good 
as  the  least  amount  of  damping  in  the  system.  Asymmetric  damping  is  not  beneficial. 

Much  like  perfect  balance,  perfect  symmetry  in  a  real  system  is  difficult  if  not 
impossible  to  achieve.  The  analysis  shows  that  small  amounts  of  stiffness  asymmetry  can 
substantially  increase  the  maximum  amplitude  of  the  rotor  response,  yet,  as  the  asymmetry 
increases,  the  amplitudes  reach  a  point  where  they  start  to  decrease  and  can  be  reduced 
roughly  to  the  amplitudes  found  with  perfect  symmetry.  If  nearly  perfect  stiffness 
symmetry  cannot  be  achieved,  it  may  be  a  better  choice  to  design  further  stiffness 
asymmetry  into  the  system  to  bring  down  the  amplitude  of  the  maximum  response.  The 
decrease  in  maximum  response  is  however  accompanied  by  an  increase  in  the  range  of 
rpms  in  which  the  amplitudes  of  the  response  is  significant. 

Cross-stiffness  terms  have  a  notable  effect  on  the  response  of  the  non-accelerating 
rotor  but  not  on  the  accelerating  rotor  unless  the  magnitude  of  the  cross  term  is  extremely 
high.  The  magnitudes  at  which  there  was  found  to  be  an  effect  in  the  model  are  probably 
unrealistically  high  for  an  actual  cross  term.  The  cross-stiffness  terms  have  very  little  effect 
on  the  amplitude  of  the  lateral  bending  response  of  an  accelerating  rotor. 

The  initial  findings  of  acceleration  scheduling  are  encouraging.  Although  the 
evaluation  performed  on  the  model  was  limited  in  its  flexibility  and  optimization,  it  did 
disclose  the  possibility  of  reducing  the  maximum  rotor  response  by  taking  advantage  of  the 
apparent  shift  in  the  critical  speed.  This  may  prove  to  be  an  extremely  promising  approach 
in  which  to  reduce  the  lateral  bending  response. 

In  this  analysis,  the  effects  of  each  of  the  parameters  were  studied  individually.  No 
analysis  was  performed  to  determine  which  combinations  of  parameters  would  give  the 
best  response.  Certainly,  a  combination  of  these  beneficial  factors  can  be  applied  to  give 
the  best  overall  reduction  in  the  maximum  lateral  bending  amplitude  of  a  rotor  transitioning 
through  its  critical  speed. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  AND  CONCLUSIONS 

An  analytical  model  was  developed  to  predict  the  response  of  a  simple  accelerating 
rotor  and  a  parametric  study  was  conducted  to  determine  possible  ways  to  minimize  the 
magnitude  of  the  response  and/or  reduce  the  number  of  cycles  occurring  at  the  higher 
magnitudes. 

Significant  reductions  in  the  maximum  response  amplitude  of  a  rotor  transitioning 
through  its  critical  speed  are  achievable.  This  is  best  accomplished  by  the  following 
methods: 

1 .  Decrease  the  forcing  function  by  balancing  the  rotor. 

2.  If  using  a  constant  rate  of  acceleration,  use  as  great  as  an  acceleration  rate  as 
practical. 

3.  Increase  damping. 

4.  Achieve  stiffness  symmetry  if  possible,  but  if  not,  use  a  great  enough 
asymmetry  to  minimize  the  response. 

Furthermore,  the  analytical  model  predicted  that  acceleration  scheduling  could 
reduce  the  magnitude  and  duration  of  the  response. 

Experimental  data  with  varying  rates  of  acceleration  were  compared  to  the  model 
with  very  accurate  results. 

B.  RECOMMENDATIONS 

There  are  many  ways  in  which  the  experimental  facility  could  be  further  modified  to 
allow  a  more  wide  variety  of  investigations  into  the  rotordynamic  field.  In  addition,  there 
are  many  ways  in  which  the  facility  can  be  further  utilized  in  it's  present  configuration. 
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1.  Determine  the  exact  time  that  each  data  acquisition  occurs. 

As  stated  earlier,  an  attempt  was  made  to  allow  for  the  exact  time  of,  or  time 
between  acquisitions  to  be  outputted  in  a  matrix  with  the  x  and  y  displacement  data.  This  is 
very  important  in  that  it  would  allow  for  phase  locking  the  data  and  being  able  to  determine 
the  precise  rate  of  acceleration  and  the  rpm  at  which  each  acquisition  occurred.  This  would 
give  a  much  more  precise  picture  of  all  aspects  of  the  rotordynamic  response.  Continued 
attempts  in  this  area  are  highly  recommended. 

2.  Increase  the  allowable  acceleration  rates  in  the  rotor  kit. 

The  rotor  kit  is  only  capable  of  a  maximum  acceleration  rate  of  15000  rpm/min. 
This  is  still  a  relatively  slow  rate  of  acceleration  and  therefore  gives  a  limited  range  in  which 
to  determine  the  effects  of  acceleration.  Discussions  with  the  Bentley  Nevada  Corporation 
revealed  the  prospect  of  changing  the  motor  controller  feedback  circuit  to  increase  these 
rates.  The  drive  motor  may  then  be  capable  of  accelerations  of  30,000  ipm/min  or  greater. 

3.  Develop  a  way  to  vary  the  stiffness  and  damping  on  the  rotor 
kit. 

The  experimental  facility  has  much  less  flexibility  than  the  computer  models  in  that 
you  cannot  vary  the  stiffness  or  damping.  As  discussed  in  Chapter  IV,  the  mounting  table 
had  a  great  effect  on  these  values.  If  the  characteristics  of  the  mounting  table  could  be 
changed,  by  shorting  the  resilient  mounts  mechanically  for  instance,  the  effects  of  varying 
stiffness  and  damping  could,  at  a  minimum,  be  qualitatively  observed. 

4.  Conduct  a  more  detailed  analysis  of  acceleration  scheduling. 

The  results  obtained  by  instantaneously  changing  the  rate  of  acceleration  discussed 
in  Chapter  III  look  promising.  Further  research  in  the  area  of  optimizing  an  acceleration 
schedule  through  critical  could  prove  beneficial.  If  analysis  does  show  that  acceleration 
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scheduling  is  a  feasible  means  of  reducing  the  maximum  amplitude  and/or  the  amount  of 
time  the  rotor  experiences  significantly  large  vibrational  amplitudes,  then  experiments 
should  be  performed  to  verify  the  analysis.  This  would  require  a  different  motor  controller 
than  is  currently  installed  on  the  experimental  facility.  A  computer  controlled  motor 
controller  capable  of  scheduling  the  acceleration  rate  would  be  necessary. 
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APPENDIX.  LIST  OF  DIMENSIONS  AND  PROPERTIES 


Rotor  Shaft: 

Material: 

Modulus  of  Elasticity: 
Length: 

Diameter 

Mass 

Density 

Rotor  Disk: 

Material: 

Width: 

Inner  Diameter 
Outer  Diameter 
Mass 


4140  Low  Alloy  Steel 

200  GPa 

45.7  cm 

0.9525  cm 

0.253  kg 

7766kg/m3 

316  Stainless  Steel 
2.54  cm 
0.9525  cm 
1.18  cm 
0.816  kg 


29.2(106)  psi 
18.0  in 
0.375  in 
0.558  Ibm 
0.281  lbm/in3 


1.0  in 
0.375  in 
3.0  in 
1.8  Ibm 
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